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Abstract 

Recent advancement in fabrication technologies enable the construction of nano-objects with 
rather rich internal structures such as double or triple quantum dots, which can then be regarded as 
artificial molecules. The main new ingredient in the study of the Kondo effect in such artificial (and 
also in natural) molecules is the internal symmetry of the nano-object, which proves to play a crucial 
role in the construction of the effective exchange Hamiltonian. This internal symmetry combines 
continuous spin symmetry SU(2) and discrete point symmetry (such as mirror reflections for double 
dots or discrete C% v rotation for equilateral triangular dots. When these artificial molecules are 
attached to metallic leads, the set of dot operators appearing in the effective exchange Hamiltonian 
generate a group which is refereed to as the dynamical symmetry group of the system dot-leads 
[mostly SO(n) or SU(n)], and the pertinent group parameters (such as the value of n) can be 
controlled by experiment. In this short review we clarify and expand these concepts and discuss 
some specific examples. In particular we concentrate on the difference between the chain geometry 
and the ring geometry. When a perpendicular magnetic field is applied in the ring geometry, its 
gauge symmetry £7(1) is involved in the interplay with the spin and orbital dynamics of the dot. 
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INTRODUCTION: KONDO MAPPING AND BEYOND 



There are numerous models in the literature of condensed matter theory, whose signif- 
icance for achieving progress in our understanding of Nature goes far beyond the original 
aim of explaining specific experimental observations. One may mention in this context the 
Bardeen-Cooper-Schrieffer's explanation of the nature of electron pairing in superconduc- 
tors, the Ginzburg-Landau equation intended for describing critical fluctuations, the concept 
of self-localization of excitations in a perfect crystal formulated by Deigen, Pekar and Toy- 
ozawa and various other seminal ideas. The explanation offered by J. Kondo for the puzzling 
shallow minimum in the temperature dependent resistivity of metals doped by magnetic im- 
purities [l| is one of the most salient examples of this kind of this scenarios. To explain it, 
consider first Kondo original idea, which was formulated within the framework of a well es- 
tablished Hamiltonian describing exchange interaction between an impurity spin S r located 
on a given site r and the spin density s r pertaining to a Fermi sea of conduction electrons at 
this site. The latter is defined by the Fourier transform of the itinerant spin Skk' = cj co .fc k / cr / 
projected on the impurity site r, namely s r = Xlkk' s kk' ex P i(k — k') • r. Here f is the vector 
of Pauli matrices for a spin 1/2. The so called sc?-exchange Hamiltonian is, 

H sd = £ kcl a c- ka + JS r ■ s r (1) 

k,cr 

where is the energy dispersion of the itinerant electrons and J is the exchange coupling 
constant. At first glance, it looks deceptively simple. However, Kondo noticed that the first 
correction to the impurity scattering amplitude of the conduction electrons beyond the Born 
approximation suffers an infrared logarithmic divergence in energy or temperature. This 
results in a singular behavior of amplitudes for an antiferromagnetic sign of the exchange 
coupling (that is, J > 0), and renders perturbation theory inapplicable below a certain 
energy scale known as the Kondo temperature. Nearly two decades of incessant efforts to 
take this singularity properly into account and to find the ground state of the system had 
crowned with finding the exact solution both numerically (in a framework of Numerical 
Renormalization Group, NRG) and analytically (by using the Bethe ansatz) j^, 4]. 

Soon after Kondo's breakthrough, theoreticians started to extend this promising con- 
ceptual framework for other physical situations and for more complex objects than simple 
localized moments. It was recognized that the Kondo mechanism should work also in sys- 
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terns exhibiting electron tunneling, where two metallic slabs are separated by thin dielectric 
layer, which forms a tunnel barrier for electrons moving from one slab to the other. It was 
shown that the magnetic impurity located somewhere near the tunneling layer plays 

the same role in tunnel conductance as magnetic impurity immersed in a metal (and subject 
to exchange interaction with Fermi sea electrons) does for impurity resistance. It was also 
shown [?| that the Friedel- Anderson model [s[ for resonance scattering of conduction elec- 
trons by the electrons occupying the 3d levels of transition metal impurities can be mapped 
on the exchange Hamiltonian (pQ) provided the strong Coulomb interaction in the 3D shell 
suppresses charge fluctuations on the impurity site. In the next stage of development orbital 
degrees of freedom were incorporated in the Kondo physics. The idea of this generalization is 
based on the fact that the magnetic impurity being put in the center of coordinates imposes 
its point symmetry on the otherwise translationally invariant crystal, and the appropriate 
description of scattered waves should exploit the formalism of partial wave expansion (ei- 
ther in spherical waves s| or in cubic harmonics [ill). Based on this idea, the generalized 
Schrieffer- Wolff model was proposed III, Il2j, where the magnetic impurity is described as an 
effective iV-component moment, but the exchange scattering is not restricted by the usual 
spin selection rule Am = 0, ±1 for the projection m of this moment. Another version of this 
model allots the impurity and band states both by spin and orbital (referred to as " color" in 
the general case) index. In case of spin s — 1/2 and N colors the symmetry of the impurity 
is SU(2N) 13|, ll4j. Appearance of additional degrees of freedom in the Kondo Hamiltonian 
might lead to the interesting scenario of overscreened Kondo effect, which arises when the 
number of conduction electron "colors" exceed that of impurity moment [ijj]. This effect 
is characterized by the non-Fermi-liquid low-temperature thermodynamics unlike the stan- 
dard Kondo effect, which only modifies (although radically)) the Fermi liquid properties of 
undoped metal [161 ]. 

Another direction of expanding the Kondo physics is realized in mapping the Kondo 
or Anderson Hamiltonian on those of other (not necessarily magnetic) systems. A great 
variety of such generalizations is possible because some quantum systems may be described 
by a pseudospin, provided their low-energy states are characterized as an effective two- 
level system (TLS) and external perturbations allow transitions between these levels. The 
first example of such generalization was suggested by the Hamiltonian describing tunneling 
between a Fermi sea electrons and an atom sitting in a double- well potential 121] - Another 
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possibility of this sort arises when the crystal field splitting is involved in formation of 
the low-energy states of impurity atom (quadrupolar Kondo effect) Q- This type of 
"exotic" Kondo system was surveyed in a detailed review [20|. One should mention also the 
possibility of involving orbital degrees of freedom in the formation of Kondo-resonance for 
an adsorbed transition metal impurity, where the orbital degeneracy is lifted by the surface 



effects, e.g. the potential of atomic step edge 



2l|. 



A powerful incentive for further extension of the realm of Kondo physics has been offered 
in 1988, when the idea of underbarrier tunneling in presence of Kondo center was extended 
on the tunneling between metallic electrodes and nanoobjects like quantum dots or small 



metallic grains 



22 



231 ] . Such nanoobject may serve as a Kondo center (a localized moment) 



provided (i) the electron spectrum is discrete due to spatial quantization, so that the level 
spacing 5e exceeds the tunneling rate, (ii) Coulomb blockade prevents charge fluctuations and 
(iii) the electron occupation number is odd, so that the effective spin of the nanoobject is 1/2. 
In this case the tunneling Hamiltonian may be mapped on the effective spin Hamiltonian (CQ) , 
and the Kondo-like singularity arises as a zero-bias anomaly (ZBA) in tunnel transparency. 
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26] in the experiments 



This theoretical prediction was confirmed ten years later [24, 
on planar quantum dots. Many experimental and theoretical studies then followed this 
experimental breakthrough ever since. 

In the course of developing this new realm of condensed matter physics, further possibil- 
ities of extending the Kondo physics were subsequently discovered and exploited. Original 
idea of using the charge fluctuations as a source of Kondo tunneling was proposed in the 
interim between the theoretical prediction of Kondo tunneling and its experimental verifica- 
tion. It was shown in this paper that in a situation, where two charge states of the quantum 
dot with occupation N, N + 1 are nearly degenerate, this dot behaves as a two-level system, 
where the fluctuating charge configuration plays part of pseudospin, whereas the real spin 
projections may be treated as channel indices. Later on this idea was further developed and 



modified 
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29]- 



Another facet of Kondo physics in nanoobjects was unveiled, when the possibility of 



Kondo effect in quantum dots wit 
several theoretical publications [30, 



i even electron occupation number was considered in 
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331 ] . In this case the quantum dot with a singlet 



ground state may become magnetically active due to externa 



either at finite energy 30[ or at finite magnetic field 31 



32 



forces, and Kondo effect arises 



331 ] . Later on it was recognized 
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that in many cases the direct mapping of the original Kondo model into such system is 
impossible, because the effective symmetry of the pertinent nanoobject is neither SU(2) 
nor SU(2N). The aim of this review is to describe various physical situations where the 
underlying nanoobject possesses complex (and in some sense unusual) symmetries which 
are characterized by non-compact Lie groups or combinations of such groups with discrete 
groups of finite rotations. 



SURPLUS SYMMETRIES 



Among the sources of surplus symmetries which enrich Kondo physics of nanoobjects one 
may find both discrete and continuous rotations stemming from complicated geometrical 
configurations of complex quantum dots, as well as those induced by external fields used in 
practical realizations of nanodevices. In this short review we will refrain from description of 



great variety of these devices, which may be found in current literature (see e.g., |34j, l35l. l36| ) . 
Fortunately, most of the relevant physics may be exposed in a relatively simple situation 
of electron tunneling through multivalley quantum dots in contact with metallic electrodes 
(leads). 

A multivalley dot is an island with electrons confined by electrostatic potential in such 
a way that the spatially quantized electrons are distributed between several valleys. These 
valleys are coupled with each other by capacitive interaction and tunneling channels. Up 
to now there are several realizations of quantum dots with two and three valleys (double 
quantum dot, DQD, and triple quantum dot, TOP, respectively). Experimentally, the first 
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391 ] . It was pointed out that these 



such realizations of DQDs go back to mid 90-es 
objects can be treated as some forms of artificial molecules with each constituent dot playing 
the role of an artificial atom 



40 



41 
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43]. 



Compared with DQD, fabrication of TQD is a much more difficult experimental task 
and the first experimental realizations of these nanoobjects appeared only recently. One 
may mention in this connection the realization of TQD with an "open" central valley 44]. 
The term "open" here means that the tunneling between the side dots and the adjacent 
leads is limited by strong Coulomb blockade (see below), whereas the central dot freely 
donates and accepts electrons to and from its own reservoir, so its role in the device is only 
to mediate indirect exchange between the two side dots. Another successful attempt to 
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fabricate a TQD was recorded in response to a theoretical proposal [45j to use this geom etry 



for realization of ratchet effect in tunneling through nanoobjects. In this realization [46] 
the charge fluctuations were suppressed by the Coulomb blockade mechanism in all three 
valleys. The feasibility of filling the TQD with 1, 2 and 3 electrons by changing the gate 
voltages was demonstrated quite recently [471 ] . 

Theoretical studies of electronics in the TQD geometry were also inspired by possible 



applications in the field of quantum information [48|, |49|. Investigation of the Kondo physics 
in TQD [5Cj was motivated by the experimental observation of molecular trimers Cr 3 on 
gold sublayers 51] by means of the tunnel electron spectroscopy, which allows to observe 



Kondo-type ZBA in conductance. Later on, other properties of these trimers such as the 



two-channel Kondo effect 



52. 



531 ] and the interplay between the Kondo effect and Aharonov- 



Bohm effect in tunnel spectroscopy 



541 ] were considered theoretically. 



Our main focus of interest is a theoretical modeling of a device consisting of a multivalley 
quantum dot, metallic electrodes and corresponding gates. The latter regulate the electron 
occupation of any particular valley and of the dot as a whole. These properties and others are 
predetermined by the geometrical position of the valleys relative to the metallic leads (source 
and drain) in the device. Basically, there are three types of such devices which are possible 
for describing electric circuits with a DQD connected with a source s and drain d. They 
may be refereed as sequential, parallel and T-shape connections (Fig. [TJi,b,c, respectively). 
The two small dots which combine to form the DQD may either be identical or may differ 
in their size. Besides, different gate voltages may be applied to different valleys. 





(a) 
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FIG. 1: Double quantum dot in sequential (a), parallel (b) and T-shape (c) geometries. Filling 
black marks the valley detached from the leads. 



In the most symmetrical configuration of identical dots the only additional symmetry 
which characterizes a DQD is the permutation symmetry P%. In close analogy with Kondo 
effect in magnetically doped metals, one should consider the total symmetry of a device 
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'source + DQD + drain'. Then the symmetry group of the sequential configuration (Fig. 
la) contains the only discrete element, namely the s-d reflection axis, and this element adds 
nothing to the P2 symmetry of the isolated DQD. The same statement is valid for the T- 
shape geometry (Fig. [Ib). In the parallel geometry the system as a whole possesses two 
reflection axes, namely source-drain s-d and left-right l-r, where the indices I and r label 
two valleys of DQD. 

Triple quantum dots (TQD) present theoreticians (and experimentalists as well) with a 
richer variety of geometrical configurations and possess more symmetries (Fig. [2]). Similarly 
to DQD these trimers may be oriented both in sequential (vertical) and parallel (lateral) 
geometries (Fig. [2h and [2b, respectively). The natural generalization of the T-shape geom- 
etry presented in Fig. lc is the cross geometry (Fig. [2b). Besides, TQD may be organized 
in a form of a triangle, which may form a closed or open element in an electronic circuit 
(Figs. [2H,e and [2]F, respectively). In the two latter cases one deals with a three-terminal 
tunnel device. We will call the conformations shown in Figs. [2fe and [2]F "ring" and "fork" 
configurations. 
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FIG. 2: Triple quantum dot in sequential (a), parallel (b), cross-shape (c), two-terminal triangular 
(d), three-terminal triangular (e) and fork (f) geometries. 

Now let us discuss the discrete symmetry elements characteristic for the above TQD 
configurations. If all three valleys are equivalent, the discrete symmetry of an isolated linear 
TQD is that of the permutation group P3. The contact with the leads adds one more 
symmetry element, the s-d reflection, provided all three dots are coupled with the leads. 
We however consider the devices, where the central dot (filled black) is not coupled directly 
with the leads. Then the P3 symmetry is lost. The only discrete symmetry element, namely 
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s-d reflection is left in a vertical geometry (Fig. Wfa), whereas both s-d and l-r reflections 
characterize a symmetry of lateral TQD (Fig. [2b). The same is valid for the cross geometry 
of Fig. [2b. The P3 symmetry is inherent in the three-terminal configuration of Fig. [2fe. 
In this case it is better to use the classification of discrete rotation group C$ v , which is 



isomorphic to the permutation group P3. One may say 5J, [55J that in the perfect triangular 



configuration the TQD imposes its C% v symmetry on the device as a whole in close analogy 



with the Coqblin-Schrieffer-Cornut version 
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121 ] of the conventional Kondo problem. In 



the geometries of Fig. [2ji,f the only element which is remained of the original symmetry 
of triangle is the l-r symmetry like in the cross geometry. External magnetic field applied 
perpendicularly to the plane of the triangle lowers the symmetry of a device by adding one 
more element, that is, chirality of the electron tunneling from the source to the drain. 

Before considering the consequences of these symmetries for the Kondo physics, several 
introductory remarks about the derivation of the Kondo Hamiltonian are in order. In an 
early period of the theoretical studies, the problem of electron tunneling through short chains 
of quantum dots under strong Coulomb blockade restrictions was formulated in terms of the 
Mott-Hubbard picture 56j, |57[]. The theory was based on the idea that electrons injected 
from the source do not lose coherence when propagating through the sequence of quantum 
dots until they leave the chain for the drain electrode. This approach is valid only for short 
enough "Hubbard chain", where the tunneling W between the adjacent valleys exceeds the 
tunneling V between the dot and the metallic leads. Using generalized Landauer method, 
one describes the tunnel transparency in terms of the Green functions of a nanoobject in 



contact with the leads [58j . Such a procedure starts with diagonalization of the Hamiltonian 
of nanoobject with subsequent calculation of renormalization of the spectrum of quantum 
dot due to tunnel contact with the leads. Early studies of these problems concentrated on 
the calculations of the Coulomb blockade peaks which arise with changing the occupation 
of the valleys (so called Coulomb staircases}. In terms of the Hubbard model, the Coulomb 
resonances are the Hubbard "minibands" [56], which arise as a result of collective Coulomb 
blockade 59] (Hubbard repulsion). It is known, however, that the spectral function of the 
Hubbard model contains also the central peak of predominantly spin origin. This peak is 
responsible for zero-bias anomalies in the tunnel conductance, which are at the center of our 
attention. 

If the "chain" contains a single dot, one deals with the conventional mapping of the 
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Anderson-like tunneling problem onto the Kondo-like cotunneling problem [22l . |23| , so that 
the central peak is indeed the famous Abrikosov-Suhl resonance pinned to the Fermi level of 
the band electrons 0, Q, 0], which is responsible for the ZBA in tunnel conductance. DQD 
in sequential geometry is the first non-trivial generalization of the single dot case, and one 
of the primary tasks is to look whether some qualitative differences from the conventional 
Kondo effect arise because of combining the features of Anderson and Hubbard models in 
the effective tunneling Hamiltonian. 



In any case, working in the above paradigm [56|, |57|], one should derive the effective 
exchange Hamiltonian H ex in accordance with the following procedure. The starting Hamil- 
tonian is chosen in the same form as the original Anderson impurity Hamiltonian [9J: 

H = H hand + H dot + H t (2) 

where i^dot describes the properties of a chain detached from the lead in terms of its eigen- 
states |A): 

A 

Here N is the operator of total electron number in the dot. The last term in ([3]) describes 
the Coulomb blockade mechanism: the total occupation of the dot M in neutral and charged 
states is fixed by the Coulomb blockade parameter Q entering the capacitive energy of a 
complex quantum dot as a whole. The eigenvalues E\ are found at fixed occupation M of 
DQD or TQD with all tunneling matrix elements W and interdot capacitive interactions 
Q' taken into account. The tunneling Hamiltonian H t intermixes the states from adjacent 
charge sectors Af, M' due to injection of extraction of an electron from the complex dot: 

^ = EEE^U aa '+h.c. (4) 

AA' b=s,d ka 

where X AA ' = |A)(A'| are the universal Hubbard operators 61]. The tunnel parameters 
are usually approximated by a single parameter V. In the Hamiltonian (J3|) these con- 
figuration changing operators describe transitions between the states belonging to different 
charge sectors (one of this sectors corresponds to the neutral CQD and another belong to 
positively or negatively charged CQD. The index b enumerates the leads (b = s,d in the 
two-terminal configurations) 
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The Hamiltonian f/band has the standard form 

#band = ^£kc[ kr7 C bka . (5) 
bkcr 

Comparing to the corresponding term in (pQ), this Hamiltonian contains one more index b. 
In the geometries with s-d reflection symmetry one may rotate the frame in such a way that 
the band operators are classified as even and odd operators relative to this reflection 

Cefccr 2 / {Csku Qifccr) ? Coka 2 ^ (c s fc,j Cd ka ) (6) 

In case of single and double quantum dots this rotation usually excludes the odd combination 



from the tunneling Hamiltonian [22|] . However in case of TQD it is not necessarily the case. 
We will return to this question in section IV. 

The cotunneling (exchange) Hamiltonian is usually obtained from ([2]) by means of the 
Schrieffer- Wolff (SW) canonical transformation [7[, which excludes the states |A) belonging 
to the charge sectors M ± 1 from the effective Fock space. At fixed M we are left solely 
with spin degrees of freedom. In a conventional situation, the relevant symmetry is SU(2) 
and the SW procedure results in an effective Hamiltonian (CD). In case of even M and 
geometries including discrete rotations, the situation is more complicated, and the SW 
procedure intermixes the states |A) belonging to different irreducible representations of the 
Hamiltonian ([3]). The corresponding terms in the effective Hamiltonians may be expressed 
by means of the corresponding Hubbard operators X AA ' . In many cases, combinations of 
these operators form closed algebras which generate non-compact groups SO(N) or SU(N) 
with N > 2, describing the dynamical symmetry of complex quantum dots. Involvement 
of these dynamical symmetries turns the procedure of mapping the tunnel problem onto 
an effective exchange problem to be more complicated than in the simpler situations which 



were briefly described in Section I. New features of the Kondo effect arising as a resu 
this procedure were described for the first time using the T-shaped DQD as an example 



t of 



43]. 



Various manifestations of dynamical symmetries in physical problems are described in the 



monographs 62j, [63] . Some mathematical aspects of the d yna mical symmetries as applied 



to the Kondo problem may be found in the recent reviews 64. 165]. 
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KONDO PHYSICS FOR SHORT CHAINS 



Short chains represented in Figs. [Th.-c, [2k-c are the most elementary objects, where many 
aspects of Kondo tunneling beyond the original paradigm H Q] of Kondo mapping may 
be demonstrated. There is much in common between the linear DQD and TQD in vertical 
and lateral geometries, although there are some effects specific only form T-shaped and 
cross-shaped configurations (Figs, [lb and [2b), which will also be emphasized in this section. 



Double quantum dots 

We start a more detailed discussion of Kondo tunneling through artificial molecules with 
the case of DQD in the vertical geometry (Fig. [Ik). Historically, this is the first general- 
ization of a single dot problem (see above). At present, the tunneling through this simplest 

reoretical discussion of Kondo physics in 



DQD is well understood. Followin g an extensive t 
vertical QD geometry (66l . 

0, HQS, El HQ 



it turned out that the most general 



description of Kondo tunneling through vertical DQD may be given in terms of 577(4) and 
5*0(4) symmetries of a low-energy multiplets in cases of odd and even electron occupation 
Af, respectively. 

As was mentioned above, the microscopic description of Kondo tunneling is analyzed in 
a framework of the generalized Anderson Hamiltonian ([21). All generic features of Kondo 
mapping are seen already in the most elementary cases of Af = 1,2. DQD in these charge 
sectors can be treated as an artificial analog of the molecular ion H^ and the neutral molecule 
H 2 , respectively. 

In case of Af = 1 the eigenstates of the Hamiltonian i7dot © are 

E 1>2 = e d T W , (7) 

where is the discrete level position in the isolated valley of the DQD and W is the inter- 
valley tunneling. These two levels forming the "Hubbard miniband" in the Af = 1 charge 
sector correspond to even and odd combinations of the electron wave functions in the double- 
well confinement potential of the DQD. We are interested here not in the resonance tunneling 
through these levels, but in the ZBA connected with the Kondo effect. The characteristic 
energy Tk which scales the Kondo effect, should be compared with the level splitting W. 
However Tk itself depends on the level splitting, so that the comparison procedure should be 
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performed self-consistently. In the limiting case Tk(W) <C W, one may ignore the odd state 
E 2 . Since the odd combination of the lead electron states in (jSJ) is also excluded from the 
problem, we immediately see that in this case the mapping procedure (SW transformation) 
reduces the tunneling problem to the case of single quantum dot with the Hamiltonian ([T]), 
where the exchange constant J is estimated as J = V 2 /Ec, and Eq is expressed via addition 
and extraction energies, i.e. the energy costs to add or subtract an electron on/from the 
quantum dot: 

Ec 1 = E- 1 + EZ 1 (8) 
E + = e d + Q ~ e F , E~ = e F - e d . 

The Kondo temperature is given by the standard equation Tk = Dexp(— l/2pj), where D 
is the characteristic energy scale for the band electrons in the leads and p is the density of 
states at the Fermi level €p. 

In the opposite limit Tk(W) ^> W, one may neglect W when calculating Tk, so that 
the dot level acquires "orbital" degeneracy. Due to this quasi degeneracy the effective 
exchange acquires an additional factor 2, so that the Kondo temperature in this case is Tk = 
D exp(— l/4pj). If fact, the difference between the two limiting cases is the manifestation 
of the S77(4) symmetry, which characterizes the spin state of an electron in the double 
well potential. This symmetry will be described more strictly when we will turn to the 
case of lateral DQD. One should note that here we encountered the first manifestation of 
the complicated structure of spin multiplets in complex quantum dots, namely with non- 
universality of Tk' it crucially depends not only on the parameters of the Hamiltonian but 
also on the effective symmetry of the low-lying spin states involved in Kondo cotunneling. 
Interpolation between two limiting cases may be described in terms of gradual SU(2) — > 
SU(4) crossover. 

The study of the charge sector M = 2 uncovers another important aspect of the Kondo 
mapping procedure. This sector corresponds to the half-filled Hubbard chain, where the 
single-electron tunneling is suppressed by Coulomb blockade (th e interior of the Coulomb 
diamond in terms of the theory of single-electron tunneling 34j, |35j). On the surface of 
it, the Kondo tunneling is also not achievable because the ground state of the DQD with 
M = 2 is a spin singlet in close analogy with the case of hydrogen molecule H 2 - However, 
the dynamical symmetry of DQD plays its part in this case as well. 
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Indeed, in case of strong Coulomb blockade (3 = W/Q -C 1 the spectrum of isolated DQD 
consists of two low-lying spin states Es,t and two charge transfer excitons E efi (even singlet 
and odd triplet), with 

E s = 2e - 2(3W, E T = 2e (9) 
E = 2e + Q, E e = 2e + Q + 2(3W . 

Only the low-energy Singlet/Triplet (S/T) pair is relevant for Kondo tunneling. Like in 
the case of M = 1, the triplet state is frozen provided the exchange gap A cx = 2/3W 
essentially exceeds Tk- In the opposite case A cx <C Tk the spin multiplet as a whole is 
involved in Kondo tunneling, and the Kondo- type ZBA may survive [69!]. One may describe 
this phenomenon in terms of the theory of conventional two-site Kondo effect [75]. According 
to this theory, the antiferromagnetic intersite exchange Jyi competes with the single-dot 
Kondo temperature Trto- At small J\% (small enough W in our case) each spin is screened 
independently and the Kondo-type ground state may be achieved. However, there exists a 
critical value J c , so that at J12 > J c the two spins are locked into a singlet state and the 
Kondo effect does not apply. 

The relevant dynamical symmetry is the 50(4) symmetry of S/T manifold ("spin rotator" 



43l|) . Again we postpone the derivation of the effective Hamiltonian for the Kondo tunneling 
and the discussion of its observable manifestations for the case of lateral DQD. Here one 
should note that in the charge sector M = 2 the exchange gap A plays the same role as the 
charge-transfer gap W in case of M = 1: the dependence TV (A) is determined by a gradual 
symmetry crossover. 

Now we turn to the parallel (lateral) geometry of Fig. [Tb. A new element, which arises 
in this case is the possibility of separate channel for each dot. The new features brought by 
this additional quantum number were discussed in several publications 0, Q. In case of 
two channels the conduction electrons retain additional "color" after rotation (Q. We will 
enumerate the states corresponding to this color by the same indices i,j = 1, 2 as those used 
for the dots because they describe the odd and even states relative to the l-r reflection like 
those in the states (JTJ),©. 

We start once more with the charge sector M — 1, where a single electron in DQD is 
distributed between two wells. It is well known that the electron in a double well may be 
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described by means of a pseudospin operator T with the following components" 

T z = J2{dld 2a -dld 1(7 ), T + = Y,dU la T- = J2d\j 2(7 . (10) 

a it a 

This vector, together with four spin vectors with the following components, 

S ij = 2 ( d h d n ~ d \l d ii) ' S ij = d h d jl S ij = d ll d rt C 11 ) 

form the set of 15 generators for the SU(4) group. 

To close the basis for the effective spin Hamiltonian, one has to introduce similar vector 
operators for the electrons in the leads, namely, the pseudospin operator t with components 



t Z — ^2 {^2ka C 2k'u Aka C Xk'a) > t+ ~ ^ C 2ka C lk'a * ~ ^ C lkcr C 2k'a- (12) 
kk',cr kk',cr kk',rr 

and four spin operators with components 

s ij = X] ( c lfcT c ifc'T _ c \ki c jk>\) > s tj = ^2 c lk^ c jk'i S ij = ^2 c lkl c jk'V ( 13 ) 

bb f bb' bb' 

r\jf\j r\jf\j A. A, 

The SW transformation carried out in terms of these operators results in the effective 
Hamiltonian, 

H c s = H hand + H dot + 2 2J JijSij • Sji + 2KT ■ t . (14) 

In the fully symmetric case which we describe here, all the effective exchange constants have 
the same value Jy — K = J. Then the Hamiltonian (fT4"j) may be reduced to a more compact 
and familiar form of exchange Hamiltonian in a fictitious magnetic field ~j\ 

H eff = # band + JS ■ S — B ■ S , (15) 

where 

s = \j2^^/3fp, s = - Yl c «^ c /3 • ( 16 ) 

a/3 kk> a/3 

Spin fermion operators f a and the matrix S involve 15 components, which are generators of 
the Lie algebra su(4), determined as 

{(t+ t~,t z , I) <g> a + , a", a z , I}- {10 1} 

where <j u (t u ) are the Pauli matrices in the spin (pseudospin) space and I is the unit matrix. 
The fictitious magnetic field has only one non-zero component, namely the t z ®I component, 

14 



and its magnitude is W, so that B ■ s = WT Z . Thus one concludes that in DQD with odd 
occupation the formally multichannel Kondo Hamiltonian (114j) may be mapped onto the 
standard Hamiltonian ( ITBT) in the case of complete channel degeneracy, and the contribution 
of additional (permutation) degrees of freedom is described in this degenerate case by a 
fictitious magnetic field 

In the charge sector M = 2 the permutation symmetry degenerates into trivial unit 
transformation. As to the spin degrees of freedom, one may try to described them by means 
of two spin operators s; and s r using the above mentioned analogy with the two-site Kondo 



problem 75j . However, such approach |77j should be used with some caution. One should 
take into account that in a situation, where both triplet and singlet two-electron states are 
involved in Kondo effect, these spins are non-independent because the kinematical constraint 
is imposed on the S/T manifold by the Casimir operator C ^ s, 2 + s 2 . In accordance with 



the prescriptions of the theory of dynamical symmetries 62|, |63|, |64j , one should construct 
two operators 

S = Si + s n R = si - s r (17) 
and impose on them the kinematic constraint 

C = S 2 + R 2 = 3 (18) 

Then the three components (S z , S + , S~) of the vector S describe the states within the spin 
triplet and transitions between them, whereas the three components (R Z ,R + ,R~) of the 
vector R describe transitions between the singlet S and the states with spin projections 
H — 1, 0, —1 of the triplet T. Six components of the vectors S,R form a closed algebra 

[S a , Sp) = ieapjSj, [R a ,Rp] = ie a p 1 S 1 , [R a ,Sp] = %e a p 1 R 1 . (19) 

and form a set of generators of 5*0(4) group. Here a, (3, 7 are Cartesian coordinate indices, 
and e Qj a 7 is the anti-symmetric Levi-Civita tensor. Two vector operators are orthogonal, 
S • R = 0. Under these constraints, two vectors r are rather fictitious than real spin 
operators. More detailed discussion of interconnections between two representations as well 
as the derivation of these operators by means of the Hubbard operators X AA> may be found 



in 43J . If one tunneling channel couples this "spin rotator" with the reservoir of conduction 
electrons, then the dynamical group SO (4) exhausts the spin degrees of freedom involved in 
Kondo tunneling. This scenario is realized in the T-shape geometry of Fig. [lb (see below). 
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More complicated is the mapping procedure for the two-channel Anderson Hamiltonian 
describing the DQD shown in Fig[Tb. In this case the discrete symmetry is explicitly involved 
in the cotunneling process, so that the SW transformation give the exchange part of if e g in 
the form 

# ex = 2 ^ JiAj ■ Sji + 2 JijRij ■ Sji + 2KT ■ t . (20) 

ij ij 

DQD with l-r reflection axis, where all the states are classified as even or odd states relative 
;o this reflection, is described by this Hamiltonian may be qualified as a "double spin rotator" 



781 ]. The corresponding dynamical symmetry group is Pi x SO (4) x 5*0(4) 



Before turning to the physical aspects of Kondo tunneling through the objects of SU (4) 
and SO (4) symmetries, one should mention that the DQD in lateral geometry with two 
channels is described by the same basic two-level Anderson Hamiltonian as the single planar 



quantum dot with two levels, one of which is occupied and anot 



model was considered in many publications 



31 



32 



review papers 



74 



33j, 



77 



rer is empty. The latter 



79 



80, 



8 II ] (see also the 



82|). So the physical manifestations of "variable" symmetry are common 



for both systems. Fig. [3] illustrates the variation of T# as a function of control parameter 
(the tunnel splitting W in case of J\f = 1 and the exchange gap A ex in case of Af = 2). In 
both cases the maximum value of Tr- correspond to the degeneracy points. 



T k (W)T k (0) 



T K (A)T K< 0) 





2T K (0) 



w 





2T K (0) 



FIG. 3: Dependence of the Kondo temperature on the control parameter W and A 
odd occupation (left panel) and even occupation (right panel). 



-A cx for 



This behavior may be easily understood within a perturbative (high-temperature) renor- 
malization group (RG) approach. This approach is based on the study of flow diagrams 
describing the evolution of effective coupling parameters Jain) due to reduction of the ef- 
fective energy scale D of conduction electron kinetic energy 



831 ] (label a enumerates the 
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vertices in the exchange Hamiltonian, 77 = In D is the scaling variable). The general form of 
the system of scaling equations is 

Here j a = pJ a are dimensionless coupling constants and c a b are numerical coefficients. These 
equations should be solved under the boundary conditions j(D Q ) = pJ a o, where -Do and J„o 
are the initial conduction bandwidth and the bare exchange integrals entering the Hamilto- 
nian if ex . When D ~ Do 3> (W, A ex ), one may neglect the splitting of energy terms in the 
system ( 12TI) . and all coupling constants J a evolve together. When the scale D is reduced 
down to the splitting energy W or A cx , the coupling parameters related to the upper level 
in the multiplets (J7|) or (Q stop to evolve, because the corresponding degrees of freedom 
are quenched below these energies. These parameters are K both in the Hamiltonians (|14j) 
and ([20]) . and the exchange vertices Jy related to the singlet-triplet transitions in (|20|) . In 
a degenerate model (fT5i) . Tr- is partially suppressed by the fictitious magnetic field (B). 
As a result the temperature T% obtained as a solution of these equations becomes explicit 
function of the splitting energy. Its behavior is illustrated in the left and right panels of 
Fig. [3] for the charge sectors Af = 1,2, respectively. In case of odd occupation Af = 1 the 
point W = corresponds to two equivalent dots decoupled from each other. The system 
remains quasi degenerate till W < T K (0) = D exp 4/ pj . With further increase of the level 
splitting, the contribution of the state E2 in Eq. (J7J) to the Kondo tunneling diminishes. In 
the asym pto tic regime |W|/Tk-(0) 3> 1, the evolution of Tk is described by the asymptotic 



equation [33|, [74 \ 



Tk(W) = (22) 



This asymptotic curves are shown bold Fig. 3. 



and detailed geometry of quantum dot 



74 



le exponent depends on the model parameter 



77 



78] 



Similar asymptotic curve describes also the evolution of T K for even occupation Af = 2 
with \W\ substituted for A ex . But in this case such behavior is characteristic only for 
A cx < 0, where the ground state is spin triplet. Tk is maximal in the "critical point" of 
S/T crossover. On the triplet side of S/T transition it diminishes due to gradual quenching 
of triplet-singlet excitations described by the operator R in the Hamiltonian ( 120]) . On the 
singlet side of S/T crossover, A cx > the Kondo temperature falls down steeply. The 
singlet-triplet crossover driven by some external parameter (magnetic field, gate voltage 
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etc) in DQD with even occupation is the most salient effect predicted and observed in these 
systems. 

Additional information unveiling specific features of Kondo screening in complex quantum 
dots can be extracted from the temperature and magnetic field dependence of tunnel conduc- 
tance G max at given W or A ex . Its behavior was discussed time and again 



71 



76 



73, 



81 



82] 



The most interesting is the behavior of conductance in the charge sector M = 2. 

On the triplet side of S/T crossover G max (T) increases with decreasing T. At high T ^> 
7V(|A ex |), it grows logarithmically, 



G max (T)/G ~ \n- 2 (T/T K ) 



(23) 



where Go — 2e 2 /h is the limiting value of tunnel conductance in a Kondo regime (unitarity 
limit, where the sum of all phase shifts on the Fermi levels 8^ equals 7r/2). This limiting 
value quantify the tunnel conductance at T — 0. According to the Friedel sum rule 0, 0] , 



generalized for the two-channel geometry 



79 



801 ] one has 



G(T = 0) = G sin -(n 2 - 



71 



(24) 



where fi(2,i) is the change in number of electrons under the Fermi level due to Kondo screening 
in two channels (even and odd). In terms of the phase shifts the argument of the sine 
function in the r.h.s. of Eq. ([SI]) reads This means that the demand of spin 

rotation invariance and l-r symmetry means that the phase shift in each channel reaches 
7r/4 at T = 0. Similarly to the theory of Kondo scattering If]], the deviation from the 
unitarity limit is of standard Fermi-liquid character 



G(T)/G = [l-n(T/T K ) 2 } 



(25) 



at T Tx. The monotonous interpolation 82| between two limiting temperature regimes 
( 123]) and (|25|) may be violated due to multistage Kondo effect. We will discuss this regime, 
when considering the case of triple quantum dots (see below). 

The influence of external magnetic field B on the tunnel conductance at low T <C B is 
easily tractable 82J: in accordance with the general theory of Kondo effect in presence of 
Zeeman splitting of the levels in the dot, in this regime T should be substituted for B in 
the asymptotic equations (123!) and (125|) . At low T -C T K additional information may be 



obtained by means of the NRG method 



81] . Due to the loss of spin rotation symmetry the 
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phase shifts 5^ become explicit functions of magnetic field. In a single channel case this 
dependence is scaled by T%, namely S a (B) = 5(0) + a(B /Tk)- In a two-channel regime, 
level splitting W between the dot levels enters this dependence explicitly. A simple equation 



describing this dependence was derived analytically 



80|, 



6 ta = 6(0) + a(B/T t ) + (-l)\B/Tl) + 0(B 3 ). (26) 



Here the parameters Tj and T/, which scale the field B depend on W. NRG calculations 81 1 
show great variety of magnitoconductance curves G(B, W). 

On the singlet side of S/T crossover the ground state of DQD is nonmagnetic, so that 
there is no room form Kondo-type ZBA at T = (see the right panel of Fig. [2D- However, 
at T > A cx the triplet state is still involved in perturbative scaling equations fl2~TT) . so 
that G(T) grows with decreasing T in accordance with Eq. fl23|) . The temperature T/S 
crossover occurs at T ~ A ex , and conductance starts to fall with further decrease of T 



761 ] . Similar effect should be observed in 



ending with exponentially small value of G(T = 0) 
the behavior of G as a function of source-drain voltage v^. At eva s > A ex the DQD shows 
Kondo tunneling, but with lowering bias the triplet channel is quenched, and conductance 
shows up a zero field dip instead of zero field peak. Apparently this type of crossover was 
observed 84j in a two-orbital planar dot with even occupation. 



T-shape geometry 

In a T-shape geometry one of two dots is detached from the leads (Fig. [lb). In the 
first experimental device of this type [38j, the role of the side (right) dot was to control 
single-electron tunneling through the left dot. Recently several new effects related to the 
Kondo-regime were discovered in this geometry. 



First, it was shown [43|] that the T-shape double dot with even occupation M = 2 
demonstrates the properties of spin rotator with SO (4) symmetry. The simplest form of the 
Hamiltonian (!20j) . namely 

H ex = 2 JS ■ s + 2 JR ■ s (27) 

was derived just for this model. Besides, it was shown that in the asymmetric T-shape DQD, 
where the Coulomb blockade in the right dot is sufficiently larger than in the left dot, the 



S/T crossover may occur because of the many-body logarithmic renormalization 
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860, 



which is determined by the renormalization group invariants E%_, namely 



E* A = E A (D)-n- 1 r A ln(nD/r / 



(28) 



The level crossing is possible because the inequality IV > for the tunneling rates T\ is 
realized for asymmetric DQD, so that the renormalization of the triplet level Et is stronger 
than that of the singlet state Eg. Similar effect may be achieved in a symmetric DQD by 
means of the gate voltage applied to the side dot (later on the possibility of inducing the 
S/T crossover by means of the the gate voltage was found also for two-orbital planar dots 

Another interesting possibility arises at the odd occupation M = 1, where the strong 
Coulomb blockade exists only in the side dot jsS, 89, 90, 91]. It was found that in this 
case the Kondo resonance arises on the background of otherwise non-correlated transport 
between the electrodes via the right dot. The interference between the resonance scattering 
and free propagation is known as Fano effect discovered originally in the optical absorption 



spectra of free atoms 92J. The Fano effect in atoms is described by the same Hamiltonian 
as the Anderson impurity hybridization 9], so there is no wonder that similar effect was 
found in the resonance impurity scattering in metals 93|]. Since the single electron transport 
through the quantum dot is also described by the Anderson Hamiltonian (jSJ), the Fano effect 
ubiquitous in resonance scattering was expected in tunnel conductance of quantum dots, and 
the structures of characteristic Fano-type form were indeed observed in the tunneling spectra 
of planar quantum dot [94]. In the latter case the Fano effect arises due to the interference 
between the resonance level in the quantum dot and the band continuum in the leads. The 
corresponding contribution to the conductance has the form 

(e + qf 



G ! Fano(e) — Gq 



e 2 + l 



(29) 



where e = 2(e— 6d))/T, q is the so called asymmetry parameter predetermined by the spectral 
characteristics of the lead electrons, T = irpV 2 is the tunneling rate. In simple terms, Fano 
effect is nothing but modification of the spectral density of conduction electrons due to its 
repulsion from the resonance level superimposed on the continuous spectrum. 

In a T-shape geometry the Fano effect arises due to superposition of the Abrikosov-Suhl 
resonance created by the Kondo effect in the right dot on the continuous tunneling spectra 



of the system 'source - left dot - drain'. As was noticed by Kang et al 



, this Fano-Kondo 
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effect looks as an antiresonance in conductance: instead of the standard Friedel-Langer 
formula fl24|) . one has 

G(T = 0) = G cos 2 |n r (30) 

where n r stands for the average occupation number of the right dot. Such "unitarity limit" 
for conductance means that in the T-shape geometry, the Kondo effect in the right dot 
results in complete suppression of resonance tunneling through the left dot and a dip arises 
in the tunnel conductance instead of the usual peak. As was mentioned above, the phase 
shifts in the Kondo regime are spin dependent, so the modification of tunnel spectra due 
to destructive Kondo-Fano interference is also spin- dependent, and the T-shape dot in this 



regime may work as a spin filter in external magnetic field [95 1. 

One more advantage of the TQD geometry is the possibility to approach the two-channel 
Kondo effect 



96 



97|. Such possibility arises in the geometry of Fig. [Tb, when the Coulomb 
blockade is strong in the left dot and the right dot is big enough so that the level spacing 
is less than the tunnel rate, but the Coulomb blockade is still strong and fixes the electron 
occupation number. In such conditions the way opens to overcome the main difficulty in 
realizing a physical system that materializes the two channel Kondo model. The necessary 
precondition for this regime is in creating two separate channels that equally screen the spin 



15J . In conventional setups an electron from one channel that hops on the dot may hop to 
the other channel and thus cause mixing between the channels. This mixing lead eventually 
to two"eigen channels" with one channel coupled stronger than the other one. The channel 
with the stronger coupling fully screen the spin and the other channel is decoupled, and 



we thus have again the single channel Kondo case. It was suggested [96] to overcome this 
mixing problem by using a large quantum dot as an additional channel. Then, the free leads 
form one channel [even one in accordance with our classification (JH])] and the large dot forms 
the second channel. The channels do not mix as transfer of electrons between them charges 
the large dot. As a result, the SW transformation maps the original Anderson Hamiltonian 
for a T-shape quantum dot onto the two-channel exchange Hamiltonian 

H 2c h = ^2 e 7. fcC U C 7,fc + S 7<A S " s 7 + B $ Z ( 31 ) 

7,fe 

Here the channel index 7 stands for the even lead-dot channel and the states in the large dot 
r. The tunnel conductance in the 2-channel non-Fermi liquid regime is realized at J2 = J r - It 
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demonstrates specific temperature dependence. In accordance with predictions of conformal 



field theory 



97, 



9Sj, 



(32) 



Practically, the zero temperature limit cannot be achieved because the Kondo screening in 
the large dot is quenched due to discreteness of its electron spectrum at T ~ Se. 

Unlike the single channel case, the magnetic field is a relevant parameter in the two 



channel case 



991 ] . Introducing the channel anisotropy parameter A c h = J r — J 2 <C J2, which 



describes deviation from the 2-channel fixed point (A c h = 0), the following equation for the 
magnitoconductance may be derived [97], 199 ] 

n 

G n fi_(T, B) = Go 



- + a sign( A ch ) — - b — In — 



(33) 



where = A&Tk/ Jf , and a, b are positive numerical coefficients of the order of 1 



Triple quantum dots 

In this section we consider new features of Kondo mapping which are characteristic for 
the chains consisting of three dots (Fig. [2h-c). Simple increase of the number of sites in the 
chain promises no new qualitative results, so we will discuss here the configurations where 
the central dot differs from two side dots in its size (and hereby by the magnitude of the 
Coulomb blockade parameter Q c ), but the latter dots are identical, so that the TQD retains 
its reflection symmetries. Both conceivable situations, namely Q c Q s and Q c Q s will 
be considered (the indices c, s are used in this section to indicate the physical quantities 
related to the central and side dots, respectively). In the first experimental realization of 



TQD [44j the former option was chosen, whereas the first theoretical study [78| was devoted 
to to the second possibility. 

In a TQD with "open" central dot and Q c Q s , its role in formation of Kondo tunnel- 
ing regime reduces to providing the channel for indirect RKKY-type mechanism of exchange 
between two localized spins formed in side dots. Thus, from the theoretical point of view 
this problem is in fact may be mapped onto that for a DQD (see Section III. A) with spe- 
cific exchange mechanism. The Kondo effect may be described in terms of two-site Kondo 
Hamiltonian, where the trend to interdot spin coupling competes with t he trend to individ- 



ual Kondo coupling between two side dots and the electrons in the leads [100 
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10 1| . One may 



mention in this connection the phase diagram of effective two-site Kondo effect calculated 



in a framework of t 
in a side geometry 



re m odel, where the pair of spins is coupled to the linear electron chain 



1021 ]. 



Qualitatively new features of Kondo mapping problem were found out in the opposite 
limiting case Q c 3> Q s . As follows from the general scheme of Kondo- mapping based on 
the dynamical symmetry of artificial molecule (Section III. A), the form of effective exchange 
Hamiltonian (120]) depends on the structure of spin multiplet of isolated multivalley quantum 
dot. Up to now only two possibilities were exploited: spin and orbital doublet for odd 
occupation Af = 1 and single-triplet pair for even occupation Af = 2. Stud y of TQD with 



78 



103]. 



odd and even occupations Af = 3 and Af = 4 give us new opportunities 

To demonstrate these opportunities let us consider TQD in lateral geometry. Genesis of 
spin multiplets may be understood from a general set-up illustrated by Fig. HI In case of 




FIG. 4: TQD in parallel geometry and energy levels of each dot e a — ev ga at W a = 0. 

Af = 3 three electrons are distributed over three dots in such a way that the state with doubly 
occupied central dot is suppressed by strong Coulomb blockade Q c . In accordance with the 
Young tableaux for a system with l-r permutation symmetry, the spin multiplet consists of 
two doublets with S = 1/2 having even and odd symmetry relative to this permutation and 
one quartet S = 3/2 with full orbital symmetry. In case of Af = 4 the spin manifold consists 
of two spin singlets S £;0 and two triplets T e o , both even and odd relative to l-r permutation. 
Varying the gate voltages v g i, v gr and playing with tunnel parameters V^ r and Wi >r , one may 
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break l-r symmetry (Fig. H]) and change the singlet-triplet splitting, so that the spin states 
are classified as Si >r and T^ r . The relative positions of energy levels in spin multiplets evolve 
as a function of model parameters and various types of level crossings occur (see [78] for 
detailed calculations). Similar situation arises for vertical configuration of Fig. [2h- 

In accordance with general theory of dynamical symmetries 64j , quasi degeneracy of low- 
lying states in spin multiplets within the energy scale ~ Tk generates special symmetries 
of TQD. For example, if the multiplet of low-lying states consists of two singlets and one 
triplet, the relevant dynamical symmetry is SO(5). If this multiplet is formed by two triplets 
and one singlet, the corresponding symmetry is SO(7), etc. The methods of constructing the 
generators for these groups are described in details in the reviews 64], |65j. As a result, an 
unique opportunity arises to change the value of index n characterizing the symmetry SO{n) 
of TQD by varying the gate voltages and other experimentally controllable parameters of 
a device. The phase diagram of vertical TQD with M = 4 calculated in 
in Fig. [5j This diagram shows great variety of phases with different symmetries from the 
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78| is presented 
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FIG. 5: Dynamic symmetries in TQD: Phase diagram in coordinates x = T;/r r , y = (ei — e c )/(e r 

Sr). 



most symmetric one P2 <8> SO (4) ® SO '(4) to conventional 5*0(3) phase where the ground 
state of TQD is spin triplet, or non-Kondo singlet ground state (shaded areas). Each phase 
is characterized by its own Tk, and this means that the ZBA in conductance should follow 
the change of the Kondo temperature, so that each crossover from one symmetry to another 
is accompanied by the abrupt change of conductance at given temperature. 
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At Af = 3 we meet a somewhat unexpected situation where Kondo tunneling in a quantum 
dot with odd occupation demonstrates the exchange Hamiltonian of a quantum dot with even 
occupation. The reason for this scenario is the specific structure of the wave function of TQD 
with N = 3. The corresponding wave functions are vector sums of states composed of a 
"passive" electron sitting in the central dot and singlet / triplet (S/T) two-electron states in 
the /, r dots. Then using certain Young tableaux [781 ]. one concludes that the spin dynamics 
of such TQD is represented by the spin 1 operator S corresponding to the I — r triplet, 
the corresponding R-operator R and the spin 1/2 operator s c of a passive electron in the 
central well. The latter does not enter the effective Hamiltonian H ex but influences the 
kinematic constraint via Casimir operator C = S 2 + M 2 + s 2 = ^. The dynamical symmetry 
is therefore SO(4) ® SU(2), and only the SO (A) subgroup is involved in Kondo tunneling. 
Similar situation, although for different reasons is realized in fork and cross geometries (see 
below). 



Remarkable symmetry reduction occurs in external magnetic field 78] • First example 



itting of S/T multiplet 



3JJ , so that the up spin 



of such reduction was found in a situation where the exchange sp 
(symmetry group SO (4)) is compensated by the Zeeman splitting c 
projection |T1) of triplet forms a pseudospin with singlet \S) and the symmetry reduction 
5*0(4) — > SU{2) takes place. In case of TQD with SO (5) symmetry, due to the same 
compensation the system may be left in a subspace {Tli, Si, S r }. The symmetry reduction 
in this case is SO (5) — > SU(3), and the Anderson Hamiltonian is mapped on a very specific 
anisotropic Kondo Hamiltonian involving only operators Rj, 

ij [iv 

where //, v are cartesian components of scalar product. Here the Kondo effect is described 
exclusively in terms of dynamical symmetry. 

Another non-standard manifestation of Kondo mapping for linear TQD is the possibility 
of two-channel Kondo effect in verti cal geometry of Fig. [2^ with Q c ^> Q s at M = 3 with 



preserved s-d mirror symmetry 1031 ] . The strong Coulomb blockade in central dot prevents 
direct s-d tunneling. The cotunneling is possible only because the wave functions of electrons 
centered on the side dots have small tails on the central dot. It is crucially important that the 
standard rotation ([6]) does not eliminate the odd channel from the tunneling Hamiltonian in 
TQD. In the situation, where the ground state of TQD is the spin doublet with even parity 
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I -D e )) the SW transformation for the original Anderson Hamiltonian results in anisotropic 
two-channel exchange Hamiltonian, 

H ex = J S S ■ s s + J d S ■ s d + J sd S ■ (s sd + s ds ) (35) 

Due to the presence of nondiagonal vertex J sd the incurable orbital anisotropy arises: 
the tunneling through two channels is controlled by the parameters J± = (J s + J d ) / 2 ± 



a/ (J s — Jd) 2 /^ + Jsd)- I n accordance with the theory of two-channel Kondo effect 15[ , this 
anisotropy makes the 2-channel fixed point unattainable, but due to strong Coulomb block- 
ade in central dot J sd =~ V 2 W 2 /Q c e 2 d is extremely small, so that one may approach the 
fixed point close enough, and the predecessor of 2-channel regime may be observed experi- 
mentally as a dip in conductance G as a function of the difference of gate voltages v gs — v gd 
applied to the side dots (this difference controls the degree of channel anisotropy). 

TQD in cross and fork geometry 

To complete the studies of linear artificial molecules we consider in this section the con- 
figurations shown in Fig. [2]^c,f). We have seen above that one may meet the situation, where 
the linear TQD with odd occupation and half-integer spin demonstrates the Kondo physics 
characteristic for even occupation with integer spin due to the fact that one of the electrons 
in the dot does not p articipate in tunneling. Here we will discuss two more mechanisms of 



such "disguise" [104 ]. 



One of these mechanisms is realized in cross geometry at occupation M = 3 under condi- 
tion Q s ^> Q d for Coulomb blockade parameters. In this case two side electrons are passive: 
the tunneling between source and drain occurs through the central dot. However, these 
passive electrons influence Kondo mapping because they are responsible for the parity of the 
3-electron wave function relative to the l-r reflection. Diagonalization of the low-energy spin 
states shows that it consists of three spin doublets and one spin quartet, and the low est st ate 
in this manifold is the doublet D u , which is odd relative to the mirror reflection (see [104| for 
details). Although the wave functions of two passive electrons do not enter explicitly in the 
indirect exchange integral arising due to SW transformation, this integral changes its sign 
due to odd parity of the state D u . Thus, in contrast to the standard paradigm of Kondo 
mapping, the effective exchange Hamiltonian corresponds to ferromagnetic coupling, which 
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is irrelevant to Kondo effect, and the Kondo-type ZBA does not arise in this case in spite 
of the fact that the net spin of quantum dot is 1/2. 

However, this is not the end of the story. The excited states in the spin multiplet which 
are Kondo active, influence the tunnel transparency and conductance at finite temperature 
and finite energies of incident electron due to non-trivial dynamical symmetry of TQD 
described above. The states involved in the Kondo effect are the even spin doublet D g 
and the quartet Q, so that the overall dynamical symmetry of TQD in cross geometry 
is SU(2) (g> SU{2) ® SU{2). We meet here the situation, which reminds the "two-stage" 
Kondo effect in DQD with TV = 2, on the singlet side of S/T crossover (see Fig. [3b and 
subsequent discussion). Here, however, there are three stages of Kondo screening, where 
the states Q and D g are quenched one after another with decreasing energy or temperature 
and the Kondo screening eventually stops at zero T. Besides, the hierarchy of tunneling 
rates Tq > T Dg > T Du exists in this charge sector, so the level crossing controlled by the 
parameters of TQD is possible in accordance with Eq. ( f2~8l) . This level crossing is shown 
in Fig. El Here the scaling variable is chosen in the form rj = \n(irD/TQ). The value of D 
is determined from the crossover condition D\ « E\(D\), where the renormalization ( 1281) 
changes for the SW regime with fixed charge, D is the initial value of scaling variable. 
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FIG. 6: Flow diagram for the levels E\ determined by the scaling invariant (|28|) . A = D u , D g ,Q, 
d = ttD/Tq. Energy is measured in meV units (see text for further explanations). 

Three points (In D u , In Dq, In D cr ,) on the abscissa axis correspond to three values of the 
control parameters where the crossover to the SW regime occurs for the ground states Ed u , 
Eq and the completely degenerate ground state, respectively By means of appropriate 
variation of the control parameters, the system may be transformed from a non-Kondo 
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regime with the ground state Ed u to the underscreened Kondo regime with the ground state 
Eq and spin S — 3/2. In accordance with the general theory of Kondo mapping (Section 
III. A), Tx is maximum in the point of maximum degeneracy. Evolution of Tr is shown in 
Fig. which should be compared with Fig. [3b. Unlike the case of singlet/triplet crossover 




FIG. 7: Evolution of Tr as a function of a control parameter D. 



in DQD with M = 2, here one deals with the crossover from the non-Kondo spin doublet 
D u to the Kondo spin quartet Q via the highly degenerate region of SU{2) (g> SU (2) ® SU (2) 
symmetry. 

Let us turn to the fork geometry shown in Fig. [2]F. In this geometry non-trivial Kondo 
physics arises already in the simplest case of odd occupation M = 1 in a situation with 



the l-r mirror symmetry. The fork may be considered as a "quantum pendulum" 49, 
Three one-electron eigenvalues are 



1051 



E D - "' 



e c =F 2W 2 /A, E Dn = e s 



(36) 



(A = max{|e s — e c \, \Q C — Q s \}) The eigenstates are classified as a non-bonding spin 
doublet D n (odd combination of the wave functions centered in the sites 1,2) and bond- 
ing/antibonding pair D^, a of corresponding even combination with the state centered in the 



cite 3 (see 



104j for details). The latter pair is the analog of resonant valence bonds (RVB) 



known in "natural" molecules. To describe this pendulum one should introduce the pseu- 
dospin vector T defined in Eq. (1101) and work with the Hamiltonian (TL4"1) . Like in the cross 
geometry, the level crossing effect as a function of control parameter takes place (Fig. [S]). 

Here Tk is nonzero on both sides of the crossover and its evolution is described by the 
bell-like curve similar to that in Fig. [3^ (although slightly asymmetric). However the 
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FIG. 8: Flow diagram for the levels E\ of the TQD in fork geometry. 



tunnel conductance is drastically influenced by the pendulum structure of the electron wave 
function. In the three-terminal fork geometry, one should consider separately the situations, 
where the bias voltage is applied between the leads 1 and 2 and between the leads 1 and 3. 
We define the corresponding components of tunnel conductance as G22 and G33, respectively. 
The Kondo anomaly in G22 is predetermined by the RVB pair, and the ZBA roughly follows 
the evolution of Tk through the crossover. More peculiar behavior is expected in 1-3 channel 
because the non-bonding state D na = 2~ 1 / 2 (d\ r7 — d\ a ) is detached from the lead 3. As a 
result the Kondo contribution to G33 manifests itself as a finite bias anomaly (FBA) in a 
situation where the ground state of TQD is Eo n . Tunnel conductance as a function of bias 
voltage in both channels is illustrated in Fig. [9j The dip in the curve b on the right panel 
reminds similar dip in the tunnel conductance of DQD with M = 2 on the singlet side of 
S/T crossover [zfjlfyj]. 
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FIG. 9: Left panel: Tunnel conductance in the channel '2—1' for D = D cr (a), D > D cr (b) and 
D < D cr (c). Right panel: Tunnel conductance in the channel '3-1' for D > D cr (a) and D < D cr 
(b). 



Thus, we see that some manifestations of the Kondo effect in TQD with odd occupation 
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may mimic those for DQD with even occupation due to specific influence of the mirror 
reflection on the structure of the electron wave functions in trimers. 



KONDO PHYSICS FOR SMALL RINGS 



In previous sections we discussed the Kondo effect in short "Hubbard chains" in contact 
with metallic reservoirs. Meanwhile, the Hubbard-like objects were studied also in closed 

Il07l |). Exp e riment ally Kondo effect on closed rings 
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109| , but the impact to theoretical 



501 ] was given by experimental studies 



ring geometries as well (see, e.g., 106, 
was observed on gutter-like dots in planar geometry 
investigation of Kondo effect in ring-like nanoobjects 

of Co trimers adsorbed on metallic sublayer jsi ] . It was shown that the basic symmetry of 
equilateral triangular triple quantum dot (TTQD, Fig. [2b) with odd occupation is SU(4) 
due to the interplay between the spin and orbital degrees of freedom, similar to that in 



two-orbital DQD (see Section I 



models triangular Co trimer 53 



LA). Special attention was paid to the case M = 3 which 



110| . In this charge sector the effective spin Hamiltonian 



contains not only exchange interaction between the spins in the dots and adjacent leads, 
but also the two-site Heisenberg exchange between spins in the neighboring dots. Magnetic 
frustrations in triangular geometry affect the spin state and therefore influences the Kondo- 
type ZBA in tunneling spectra. These spectra were calculated by the numerical RG and 
quantum Monte-Carlo methods. Besides, it was found 53j] that in case of complete channel 
isotropy the non-Fermi-liquid regime arises from the interplay of magnetic frustrations and 
Kondo physics. 

Another phenomenon, which interplays with the Kondo physics is the Aharonov-Bohm 
oscillation of tunnel transparency in magnetic field directed perpendicularly to the plane of 
triangle. This effect may be seen already for the Hubbard ring with 1/3 occupation (TTQD 
with M = 1), where there is no room for exchange interaction between spins localized in 
neighboring sites and concomitant magnetic frustrations. The starting point for solving; 



the problem of interplay between Kondo and Aharonov-Bohm phenomena 0, lllll is the 
Anderson Hamiltonian (j2J) rather than the exchange Hamiltonian ffl4|) . 

In accordance with the general scheme discussed in Section II, one should start with 
the diagonalization of the Hamiltonian of 3-site Hubbard ring. The point symmetry of this 
equilateral triangle is C% v . This group describes the "orbital" degrees of freedom, whereas the 
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continuous spin symmetry is usual SU{2) symmetry of spin 1/2. According to the irreducible 
representations of C% v group, the spectrum of TTQD consists of three levels A = DA, DE±. 
Here as usual D stands for spin doublet, A is the fully symmetric orbital singlet and E± are 
two components of orbital doublet. The energies of these states in out-of-plane magnetic 
field B are 

E Dr (p) = t-2W cos (p - | J . (37) 

such that for negative W and for B = 0, p = 0, 27r/3, 4n/3 correspond respectively to 
A = A, E± with the ground state DA, so that the orbital degrees of freedom are quenched at 
low temperature and the SW mapping procedure ends with conventional Kondo Hamiltonian 
([T]). However, using the magnetic field as a control parameter, one may initiate level crossing 
by varying the magnetic flux $ through TTQD. This level crossing is shown on the upper 
panel of Fig. [TUJ In each crossing points and its nearest vicinity the dynamical symmetry 
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FIG. 10: Upper panel: Evolution of the energy levels Ea (solid line) and E± (dashed and dash- 
dotted line, resp.) Lower panel: corresponding evolution of conductance (Go = ire 2 /h). 



of TTQD is SU(4), so that the symmetry crossovers SU(2) — > SU(4:) — > SU(2) occur at 
$ = (n + |)<I>o, where $o is the quantum of magnetic flux. Each crossover is accompanied by 
the change of Tk from exp(— 1/2 J) to exp(— 1/4 J) and back (see Section UFA). The ZBA 
peak in the two-terminal tunnel conductance changes accordingly (Fig. [101 lower panel). 

Although formally there are three tunneling channels, the non-Fermi liquid Kondo regime 
cannot arise, because the non-diagonal components Jy appear in the exchange Hamiltonian 
(|14p . Further diagonalization should be done by means of rotating frameworks for Bloch 
electrons. This diagonalization introduces irremovable channel anisotropy, so that the 2- 



channel non-Fermi-liquid fixed point cannot be achieved unlike the case considered in 



53|, 
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where the channel isotropy was postulated from the very beginning. 

In order to realize the Aharonov-Bohm interference, one should use the two-terminal 
geometry shown in Fig. [2U. In this case there are two paths (1-3) and (2-3) for single 
electron tunneling between source and drain. Interference of two waves in the drain results 
in famous Aharonov-Bohm oscillations. The field B affects the lead-dot hopping phases. 
In the chosen gauge, the hopping integrals are modified as, W — > W / exp(i$i/3), — > 
V r s exp[±i($i/6 + $ 2 /2)], where $1,2 are magnetic fluxes through the upper and lower loop 
of the device. As a result the exchange Hamiltonian reads 

H = J S S ■ s s + J d S ■ s d + J sd S ■ (s sd + s ds ) + KT ■ t (38) 

(the latter term becomes actual when the magnetic field induces level crossing in accordance 
with Fig. [TU]). Magnetic flux enters the coupling constants J s , J d , J sd , K via SW trans- 
formation. As a result the constant J sd (&i, $2) turns into zero at some values of magnetic 
flux, so that the Aharonov-Bohm interference completely blocks Kondo transparency. Thus 
TTQD serves simultaneously as a Kondo "pass valve" and as an Aharonov-Bohm interfer- 
ometer. It should be stressed that both the continuous spin degrees of freedom and discrete 
"rotations" of triangle are involved in these two phenomena in TTQD. 

CONCLUDING REMARKS 

Among many aspects of Kondo tunneling through complex quantum dots we have chosen 
for this review only the symmetry related properties predetermined by the structure of 
the low-lying states in the spin multiplet characterizing the fixed charge sector of complex 
quantum dot. New features, which are introduced by the dynamical symmetries in the 
Kondo physics are the multistage process of Kondo screening, symmetry crossovers driven 
by experimentally tunable control parameters, interplay between continuous spin rotation 
symmetry and discrete point symmetry of nanodevices. The main tool of experimental 
monitoring of variable symmetries is study of temperature and magnetic filed dependence 
of zero- and finite-bias anomalies in tunnel conductance. 

Among other facets of Kondo effect in nanostructures one should mention non-equilibrium 
Kondo effect at finite bias and under light illumination, where both quantum dots and leads 
are far enough from thermodynamic balance. Under this conditions such phenomena as 
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spin relaxation, dephasing and decoherence influence the tunnel transport in Kondo regime. 
Real atoms and molecules also may be included in electric circuit by means of advanced 
experimental techniques (scanning tunnel spectroscopy, break-junction method etc). In this 
case phono n- and photon assisted processes should be taken into account, which result in 
interplay of Kondo resonance tunneling with various "polaronic" and "excitonic" effects. 
Besides, Kondo processes may be included in nano-electro-mechanical shuttling, so that 
Kondo physics of movable objects becomes one of challenging items on the agenda. 
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